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Abstract 

We derive the ratio of total blue luminosity to total baryon mass, Lb /Mi,, for massive (Mgas at the Abell 
radius is > 1 x clusters of galaxies up to z ~ 1 from the literature. Twenty-two clusters in our 

sample are at z > 0.1. Assuming that the relative mix of hot gas and galaxies in clusters does not change 
during cluster evolution, we use Lb/Mi, to probe the star formation history of the galaxy population as a 
whole in clusters. We find that Lb /Mb of clusters increases with redshift from Lb /Mb — O.O24(Lb/M)0 
at 2; = to ~ O.O6(Lb/M)0 at 2: = 1, indicating a factor of 2 — 3 brightening (we assume Hq = 70 km s~^ 
Mpc"^). This amount of brightening is almost identical to the brightening of the M/Lb ratio of early-type 
galaxies in clusters at 0.02 < z < 0.83 reported by van Dokkum et al. (1998). We compare the observed 
brightening of Lb /Mb with luminosity evolution models for the galaxy population as a whole, changing 
the e-folding time of star formation r by 0.1 < r < 5 Gyr and the formation redshift zp by 2 < zp < 00. 
We find that r = 0.1 Gyr 'single burst' models with zp > 3 and r = 5 Gyr 'disk' models with arbitrary 
zp are consistent with the observed brightening, while models with r = 1 — 2 Gyr tend to predict too 
steep brightening. We also derive the ratio of blue luminosity density to baryon density for field galaxies, 
adopting flbh^ = 0.02, and find that blue luminosity per unit baryon is similar in clusters and in fields up 
to z ~ 1 within the observational uncertainties. 
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1. Introduction 

Clusters of galaxies are suitable objects for study- 
ing the evolution of galaxies in the dense environments. 
Recent observations based on large telescopes including 
Hubble Space Telescope (HST) have been revealing the 
morphology-dependent evolution of galaxies in clusters 
up to z ~ 1. The evolution of elliptical and SO galaxies 
has been found to be reproduced well by the so-called 
single burst model (e.g., Schade, Barrientos, Lopez-Cruz 
1997; Elhs et al. 1997; Kodama et al. 1998; van Dokkum 
et al. 1998). Schade et al. (1996) found that spiral 
galaxies in clusters brighten by ~ 1 mag with redshift up 
to z ~ 0.5, and that this brightening is similar to that 
of field spiral galaxies in the same redshift range (see, 
however, Vogt et al. 1997 and Lilly et al. 1998 for the 
evolution of field spiral galaxies). Morphological studies 
based on HST imaging suggest that a transition from 
spiral galaxies to SO galaxies may have occurred in clus- 
ters since z ^ 0.5 (Dressier et al. 1997). An attempt of 
measuring the star formation rate of individual galaxies 
has also started (e.g., Balogh et al. 1998; Poggianti et al. 
1999). 



In this paper, we study the global (or average) star 
formation history of the galaxy population as a whole in 
clusters of galaxies. To do so, we derive the ratio of to- 
tal blue luminosity (Lb) to total baryon mass {Mb) for 
clusters of galaxies up to z ~ 1. The quantity Lb/ Mb, 
including its evolution, should reflect when (and what 
fraction of) baryons {— primordial gas) in clusters are 
converted into stars. Similar studies have been done for 
field galaxies. The global luminosity density I [L@Mpc~^] 
in various wavelengths has been measured by many work- 
ers on the basis of observations of field galaxies (e.g., Lilly 
et al. 1996; Madau, Pozzetti, Dickinson 1998). If the 
density parameter of baryons (lib) is given, one can com- 
pute from / the mean luminosity per unit baryon mass in 
fields (l/ph)- We will derive J5-band l/ph in fields up to 
z ~ 1 and compare it with Lb /Mb of clusters. 

The structure of this paper is as follows. In section 2, 
we present the data of nearby and distant clusters used 
to derive Lb /Mb- We compare Lb /Mb with predictions 
of simple luminosity evolution models in section 3. A 
comparison with 1/ ph in fields is also given in section 3. 
We summarize our conclusions in section 4. 

We adopt h = 0.7, VLq = 0.2, and Aq = throughout 
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this paper unless otherwise stated, where h is the Hubble 
constant in units of 100 km s~^ Mpc^^, Hq is the density 
parameter, and Aq is the cosmological constant. Under 
this assumption, the present age of the universe is 11.8 
Gyr. The value h = 0.7 is taken from recent determi- 
nations of ffo (e-g-) Freedman 1999). To adopt different 
values for and Aq in the observationally reasonable 
ranges of 0.2 < fJo < 1 and < Ao < 0.8 does not 
significantly change our results. 

2. Data 

We divide clusters of galaxies into nearby clusters {z < 
0.1) and distant clusters {z > 0.1). We assume that 
evolutionary effects are negligible for z < 0.1 clusters 
and regard their properties as those of the present-day 
clusters. 

We adopt the B band to measure the luminosity of 
galaxies, and think that it is the best compromise. Clus- 
ter luminosities have been measured mainly in optical 
bandpasses such as B, V, or R. Among the optical band- 
passes, B is most sensitive to the luminosity evolution of 
galaxies. Though ultraviolet wavelengths such as U are 
much better for measuring star formation, data in such 
wavelengths are very few. As for field galaxies, there are 
a lot of measurements of luminosity density of field galax- 
ies in the B band, which enables us to compare Lb /Mb 
with iB/pb- 

2.1. Nearby Clusters 

We use the sample of nearby clusters given in Arnaud 
et al. (1992) to derive Lb/Mi, of the present-day clusters. 
Arnaud et al. (1992) compiled a sample of 27 clusters of 
galaxies, where total V^-band luminosity (iv), morpho- 
logical type mix of galaxies (E, SO, and S), and gas mass 
within a radius of 1.5h~^ Mpc (the Abell radius) are 
given. Morphological type mix is available for 18 clus- 
ters. 

Arnaud et al. (1992) found in their clusters a strong 
dependence of -LvA^gas on Ly: iy/Afgas oc L^'^'^ . If 
such a strong dependence holds in the whole mass range 
of clusters, it would very much complicate a compari- 
son of Lb /Mb among clusters having different masses. 
Thus, we first examine for what clusters such a strong 
dependence exists. Figure 1 plots Ly against Mg^s for 
all the clusters in Arnaud et al. (1992). The dependence 
found in Arnaud et al. (1992) is shown as the dashed 
line. The solid line, on the other hand, is a regression 
hue between Ly and Mgas for Mgas > 1 x lO^^/i^^-^M® 
clusters, Ly oc M^^. This is close to a linear regres- 
sion, i.e., a constant Ly/M^^, indicated as the dotted 
line. Thus, the strong dependence of Ly/^^gas on Ly (or 
equivalently on Mgas) found by Arnaud et al. (1992) is 
probably due to the inclusion of less massive clusters. 
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Figure 2 shows iy/Afgas as a function of the frac- 
tion of luminosity emitted from elliptical and SO galax- 
ies to the total luminosity for 18 clusters with type 
mix data. It is found that iy/Mgas is constant for 
Mgas > 1 X l^y^^hr"^-^ A4(^ clusters irrespective of type 
mix. This implies that Ly/Mgas is not sensitive to the 
change in the populations of galaxies for massive clusters. 

Figures 1 and 2 demonstrate that the dependence of 
Ly/Mgas on Ly (or on Mgas) is much weak for massive 
clusters. This is also supported by Renzini (1997) who 
found that rich clusters have a fairly constant Mgas to 
S-band luminosity ratio. The reason why poor clusters 
have a relatively higher Ly/Mgas value is not clear, but 
a possible explanation is that a significant fraction of hot 
gas in poor clusters has escaped from the clusters during 
cluster evolution owing to their shallow gravitational po- 
tentials, resulting in a higher Ly/Mgas value (e.g., Ren- 
zini 1997). 

In any case, in what follows we assume that Lb /Mb 
is constant for Afgas > 1 x 10-^^/i~^-^Mg, clusters and use 
them to derive the average Lb /Mb of nearby clusters. We 
will see in the next subsection that all the distant clusters 
adopted in this paper have Mgas > 1 x 10^"^/?-~^ ''^A4g). This 
promises a fair comparison of Lb /Mb between nearby and 
distant clusters. 

For Arnaud et al.'s (1992) clusters which have mor- 
phological type mix, we compute total B-h&nd luminosity 
(Lb) from Ly using B-V = 0.96 (E), 0.85 (SO), and 0.68 
(S) (see Fukugita, Shimasaku, Ichikawa 1995). For clus- 
ters without type mix data, we adopt B — V = 0.85 ± 0.2 
as the average color of galaxies. We compute baryon 
mass Mb from gas mass Mgas using: 

Mb = Mgas + (M/Lb)*Lb, (1) 

where (M/Lb)* is the mean mass-to-luminosity ratio of 
the stellar population in galaxies. We neglect atomic 
and molecular gas in galaxies. We adopt (M/Lb)* = 
(6 ± 3)/i(M/Lb)0, which roughly covers the mass-to- 
luminosity ratio of elliptical galaxies (van der Marel 1991; 
Pizzella et al. 1997) and of spiral disks (Bahcall 1984; 
Broeils, Couteau 1997). The errors in Lb /Mb contain (i) 
errors in Ly which are given in Arnaud et al. (1992), (ii) 
errors in the mean B ~V (only for clusters without type 
mix data), and (iii) errors in (M/Lb)* in equation (1). 
Since most of the baryons in clusters are in form of hot 
gas, the error in Mb due to (iii) is only about 5 %. 

Figure 3 presents Lb/ Mb as a function of Mgas- The 
filled and open circles indicate Arnaud et al.'s (1992) 
clusters with and without type mix data, respectively. 
As seen in figure 1, a clear trend is seen in figure 3 that 
clusters with Mgas ^ 1 x 10^^h~^-^Mg, have systematically 
higher Lb /Mb. To derive Lb /Mb of nearby clusters, we 
not only remove clusters with Mgas < 1 x lO^^/j-^-^M® 
but also remove clusters without type mix data because 
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the uncertainties in Lb of these clusters are on the aver- 
age larger than those for clusters having type mix (To in- 
clude the clusters without type mix data hardly changes 
the result, though). 

Twelve out of the 21 clusters with Afgas > 1 x 
lQ^^h~'^-^M^ have type mix data, and their mean Lb /Mb 
for h = 0.7 is 

Lb/Mi, = (0.024 ± O.OO4)(Lb/M)0, (2) 

which we regard as the representative value for the 
present-day clusters. The contribution from elliptical 
and SO galaxies to total B luminosity is on the aver- 
age (69 ± 13)% for the 12 clusters, implying that these 
clusters are dominated by early-type galaxies (See figure 
2). Clusters which have Mgas < 1 x 10^^/i~2-^M® tend to 
be less dominated by early-type galaxies. For example, 
the Virgo cluster, which has Mgas = 0.44 x 10^^/i~^-^M®, 
has Lb{E + S0)/LB(tot) = 45%. 

For the 12 clusters, we estimate the ratio of stellar mass 
to baryon mass to be M^/Mb = 0.10 ± 0.05 {h = 0.7), 
using (M/Lb)* = (6 ± 3)/i(M/Lb)0. This moans that 
only ~ 10 % of baryons have been used to form stars to 
date in rich clusters. 

2.2. Distant Clusters 

Searching for total luminosity and gas mass data of dis- 
tant clusters in the literature, we take 22 clusters, among 
which thirteen are from the CNOC cluster sample (Carl- 
berg et al. 1996; Lewis et al. 1999). We do not apply any 
selection criterion to compile our sample. Data of these 
clusters are given in table 1. All the clusters have rest- 
frame luminosity (either B, V, or r band) and gas mass 
measurements. For a cluster whose rest-frame luminos- 
ity is in the V ov r band, we use an observed or assumed 
color to convert the luminosity to the rest-frame B-band 
luminosity. 

When computing Mb from Mgas by Mb = Mgas + M^,, 
we use Mi, /Alb = 0.10 {h = 0.7) which is the value for the 
nearby clusters. Distant clusters may have lower M^^/Mb 
values than the nearby clusters, but the uncertainties in 
Lb /Mb due to this effect are at most ~ 10%, which is 
negligible for our discussion. Below are the references to 
the nine clusters and the CNOC clusters. 

(i) Abell 1413 [z = 0.14) and Abell 1689 {z = 0.18) 

These two clusters are taken from Cirimcle, Nesci, 
Trevese's (1997) sample. This sample consists of 12 clus- 
ters with 2; < 0.2 for which Mgas and iy within a radius 
of 0.75ft-~^ Mpc are given. We calculate Lb assuming 

rest-frame B -V = 0.85 ± 0.2. 

(ii) Abell 2218 {z = 0.18) 

We adopt Mgas and Lb from Squires et al. (1996). 

(iii) Abell 2163 [z = 0.20) 



Mgas and Ly are given in Squires et al. (1997). Lb is 
computed assuming rest-frame B — V = 0.8 ±0.2. 

(iv) CL0500-24 (z = 0.32) 

Mgas is taken from Schindler and Wambsganss (1997). 
Ly is given in Infante et al. (1994), and Lb is computed 
using rest-frame B — V = 0.85 (The mean apparent color 
of this cluster, V — I = 1.8, reported by Infante et al. 
1994 corresponds to B — y = 0.85 in the rest frame). 

(v) CL0939-I-47 (z = 0.41) 

Schindler et al. (1998) foimd two substructures in this 
cluster, implying that this cluster has not been virialized 
yet. Mgas and Lb adopted here are the sum of the values 
for the two substructures given in Schindler et al. (1998). 
Dressier et al. (1997) reported the fraction in number of 
elliptical and SO galaxies to be 55%. 

(vi) RXJ1347.5-1145 {z = 0.45) 

Mgas is taken from Sahu et al. (1998) and Lb is com- 
puted from Mtot and Mtot /Lb given in Fischer and Tyson 
(1997), where Mtot is the total mass of a cluster. 

(vii) CL0016-I-16 {z = 0.55) 

Mgas is taken from Neumann and Bohringcr (1997) and 
Lb is computed from the r-band total luminosity given 
in Carlberg et al. (1996) using rest-frame B — r = 0.97, 
which corresponds to the observed g — r color of 1.455 
(Carlberg et al. 1996). Dressier et al. (1997) reported 
the fraction in number of elliptical and SO galaxies to be 
73%. 

(viii) AXJ2019+1127 {z = 1.01) 

Mgas is taken from Hattori et al. (1997), and Lb is com- 
puted from Ly (Benitez et al. 1998) assuming rest-frame 
B — V = 0.7±0.3, which roughly covers the expected col- 
ors of elliptical and spiral galaxies at z = 1. 

CNOC Clusters 

Carlberg et al. (1996) give rest-frame r-band luminosity 

at the virial radius for 16 clusters at 0.17 < z < 0.55. 
Out of them, 14 clusters have gas mass measurements 
(figure 4 of Lewis et al. 1999). Since the maximum ra- 
dius at which gas mass is plotted, r = 600h~^ kpc for all 
the clusters, is smaller than the virial radii, we derive for 
each cluster the r-band luminosity at 600/i~^ kpc from 
the value at the virial radius assuming that luminosity is 
proportional to radius. The ratio of 600/i~^ kpc to the 
virial radii of 14 clusters is on the average 0.48, implying 
that a large factor of conversion is necessary. For each 
cluster, we then transform the rest-frame r-band lumi- 
nosity into the rest- frame B-band luminosity on the basis 
of the observed g — r color given in figure 5 of Carlberg 
et al. (1996). CL0016-M6 is among the 14 clusters. As 
seen in (vii), we have adopted for gas mass of this clus- 
ter the measurement given in Neumann and Bohringcr 
(1997) because their value is at r = 1.67h~^ Mpc, which 
is very close to the virial radius where r-band luminosity 
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is measured. The number of the CNOC clusters adopted 
here is thus 13. 

The radii used for measuring Lb /Mb differ among the 
clusters and most of them are smaller than the Abell ra- 
dius (See table 1). Unfortunately, it is not clear whether 
Lb /Mb measured at these small radii represent global 
values, i.e., values at the Abell radius, though there is 
a study that the Ls/Mf, of the Coma cluster is nearly 
constant between a radius of ~ 0.4:h~^ Mpc and the 
Abell radius (Taguchi et al. 1999). In this paper, we 
assume that the values of Lb /Mb derived here represent 
the global values of individual clusters (For CL0939+47, 
see, however, the next section). 

Figure 4 plots Lb against Mgas for the 12 nearby clus- 
ters and the 22 distant clusters. Both Lb and Mgas are 
values at r = 1.5/i~^ Mpc, which are derived from raw 
values on the assumption that Lsir) and Mgas(r) are 
proportional to radius r. The thick solid line indicates 
the best fit of a linear law. Lb oc Mgas, to the nearby 
clusters. The thin solid line and the dotted line corre- 
spond to a similar fit to distant clusters at 0.1 < z < 0.4 
and 0.4 < z < 0.7, respectively. It is found that the 
average B luminosity at a given gas mass increases with 
redshift. This should reflect some evolution of Lb /Mgas. 
Note that the range of gas mass is similar between the 
nearby and distant clusters: there is no distant cluster 
in our sample whose gas mass at r = 1.5h~^ Mpc is less 
than 1 X 10^^/i~^-^M®. Note also that no clear depen- 
dence of Lb /Mgas on Mgas is seen either in the nearby 
cluster sample or in the distant cluster sample. 

We mention here the effects of changing and Aq on 
estimates of Lb /Mb for distant clusters. Let dL{z) and 
dA{z) be the luminosity distance and the angular diam- 
eter distance to a cluster at z, respectively. The ratio 
Ls/Mgas is proportional to (^^^'^{z), because of Lb oc 
d\{z), Mgas cx d\^{z), and dL{z) = (l + z)2rf^(z). Thus, 
LB/Mgas depends on fio and Aq through d']^^'^{z). Since 
Mb is dominated by Mgas, the dependence of Lb /Mb on 
f2o and Aq is very close to that of LB/Mgas- Figure 5 
shows the ratio of LB/Mgas(rio, Aq) to LB/Mgas(fio = 
0.2, Ao = 0) as a function of redshift for two sets of 
(ilo, Ao). Since dL{z) is a decreasing function of f7o and 
an increasing function of Aq up to at least z = 1, we find 
from this figure that the change in Lb /Mgas due to the 
change in fio and Ao in the ranges of 0.2 < f^o < 1 and 
< Ao < 0.8 is less than ±10% for z < 1 clusters, which 
is negligible for our discussion below. 

3. Results and Discussion 

Figure 6 shows Lb /Mb of clusters as a function of 

redshift. The filled circles present the distant clusters 
and the filled square indicates the average Lb /Mb of 
the nearby clusters. Clusters without errors (but for 



CL0939-I-47 and CL0016-I-16) are the CNOC clusters. 

Lines indicate model predictions, which will be discussed 
in the next subsection. It is found that Lb /Mb increases 
with z, from LB/Mb = 0.024(Lb/M)q at z = to 
~ 0.06(Lb/M)o at z = 1, though the error in each dis- 
tant cluster is fairly large. CL0939-I-47 deviates largely 
from this trend. We suspect, however, that the observed 
Lb /Mb of this cluster does not represent the real, global 
value, because the observed LB/Mb is the value for two 
substructures whose radii are only r = 0.14/i~^ Mpc. The 
CNOC clusters seem to have a larger scatter in Ls/Mb. 
This may reflect uncertainties due to a large factor of 
the conversion of optical luminosity from the value at 
the virial radius to that at r = QQQh~^ kpc. 

There are two opposite explanations for the increase in 
LB/Mb with redshift. One is the brightening of Lb due 
to the luminosity evolution of the galaxy population as a 
whole. Note that galaxy mergings, even if they occur, do 
not change the total mass of the galaxy population and 
that star formation which could be triggered by mergings 
can be treated in the framework of the 'pure' luminosity 
evolution of the galaxy population. The other explana- 
tion is that the mass of baryons (~ hot gas) per galaxy 
decreases with redshift. However, this explanation seems 
to be less plausible, because no significant evolution has 
been observationally found for the global properties of 
clusters at z<l (e.g., Schindler 1999) (This result is, 
however, mainly for X-ray properties, and the evolution 
of the galaxy distribution in clusters is not well known) . 

In what follows, we take the former explanation as our 
hypothesis, i.e., we assume that the increase in LB/Mb 
found here is due to pure brightening of the galaxy pop- 
ulation as a whole. Then the increase found here corre- 
sponds to brightening by ~ 1 mag of the galaxy popu- 
lation. In the next subsection, wc compare the observed 
brightening with predictions of simple luminosity evolu- 
tion models of galaxies. 

3.1. Comparison with Luminosity Evolution Models 

We characterize the evolution of Lb, the i?-band lu- 
minosity summed over all galaxies in a cluster, by two 
parameters: the star formation timescale r and the for- 
mation redshift zp. In other words, we assume that all 
galaxies are formed at the same redshift zp and that the 
e-folding time of star formation summed over all galax- 
ies is T (Gyr) . We compute B-band luminosity using the 
population synthesis code developed by Kodama and Ari- 
moto (1997). The values of r and zp examined here are 
T = 0.1, 1, 2, 3, and 5 Gyr and zp = 2, 3, and oo. Models 
with r = 0.1 Gyr correspond to elliptical galaxies and 
models with r = 5 Gyr are for spiral disks like that of 
our Galaxy. We do not examine r > 5 Gyr, since to take 
T > 5 Gyr leads to too blue colors at z = which are in- 
consistent with the observed colors of galaxies in nearby 
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clusters (The mean B — Rc of the Virgo and Coma clus- 
ter galaxies is 1.2 and 1.8, respectively [Andreon 1996 for 
Coma and Young and Currie 1998 for Virgo], while the 
models with t — 0.1 and t — 5 Gyr give B — Rq = 1.6 
and B — Rc = 1.1, respectively, at an age of 12 Gyr). We 
also set the lower limit of zp to be 2 following the tra- 
ditional pure luminosity evolution models which assume 
that elliptical/so galaxies and spiral galaxies are formed 
at high redshifts (z > 2) and which broadly succeed in 
reproducing observed properties of these galaxies (e.g., 
Kodama et al. 1998; Shimasaku and Fukugita 1998). 

Figure 6 compares the observed Lb/M^ with predic- 
tions. Predicted values of L-g/Mi, are normalized to 
match the observed value at z = 0. In other words, mod- 
els are used to predict relative brightening (or fading) of 
Lb as a function of z. Panels (a), (b), and (c) are for 
zp — oo,3, and 2, respectively. From panel (a), we find 
that all the models reproduce the observation. If, how- 
ever, zp = 3 is adopted (panel [b]), models with t = 1 
and 2 Gyr give too steep brightening compared with the 
observation. This trend is strengthened for the zp = 2 
case (panel [c]): the r = 0.1 Gyr model also becomes in- 
consistent with the observation, though the discrepancy 
is at less than 2a levels. The allowed range for r is de- 
pendent on zp, and we cannot rule out any value of r 
on the basis of the current data if we permit zp = oo, 
though 'single burst' models with r = 0.1 Gyr and 'disk' 
models with t — 5 (and t — 3 Gyr models) match the 
observation for a wider range of zp toward lower redshifts 
than the other (r = 1, 2 Gyr) models. 

It is interesting that 'single burst' (r = 0.1 Gyr) models 
and 'disk' (r = 5 Gyr) models are consistent with the 
observation. In this paragraph, we concentrate on these 
models, and examine which are more consistent with the 
observed luminosity evolution of individual galaxies in 
clusters. Various observations suggest that elliptical and 
SO galaxies in clusters brighten by ~ 1 mag from z = 
to z = 1, which is consistent with the single burst model 
(e.g., Schade, Barrientos, Lopez-Cruz 1997; Ellis et al. 
1997; Kodama et al. 1998; van Dokkum et al. 1998). 
Spiral galaxies have also been found to brighten by ~ 1 
mag (e.g., Schade et al. 1996), though observations are 
not as many as those of elliptical and SO galaxies. The 
amounts of brightening of E/SO and spiral galaxies are 
similar to each other, and they are in agreement with 
the predictions of t = 0.1 Gyr models (with zp > 3) 
and T = 5 Gyr models. However, our distant clusters 
are rich clusters and thus it is likely that elliptical and 
SO galaxies dominate in these clusters. Hence, the r — 
0.1 Gyr models seem to be more plausible for describing 
the evolution of Lb /Mb. This is also supported by the 
fact that the mean color of galaxies in the Coma cluster, 
which is a very rich nearby cluster and is likely to be a 
counterpart of the distant clusters studied here, agrees 
with the color predicted by the r = 0.1 Gyr models. 



Van Dokkum et al. (1998) present observations of Al / L 
ratio in the B band of early-type galaxies in five clusters 
at 0.02 < z < 0.83. They find that the M/L ratio evolves 
as AlogM/Le oc -0.40z {flo = 0.3, Aq = 0), which is 
consistent with single-burst models with Zp > 1.7—2.8. If 
the evolution of M/L found by van Dokkum ct al. (1998) 
is understood as pure luminosity evolution of Lb , the for- 
mula AlogM/ie oc — 0.40z implies that Lb brightens 
by a factor of 2.5 from z = to 1, which is in excellent 
agreement with the brightening of Lb/Mjj found in this 
study. Measurements of M/L of individual galaxies are 
a direct measurement of the effects of luminosity evolu- 
tion occurred in galaxies, while measurements of Lb /Mb 
of clusters are less direct. Note, however, that informa- 
tion contained in Lb /Mb is different from that in M/L of 
individual galaxies. The quantity Lb /Mb describes the 
evolution of luminosity summed over all galaxies in clus- 
ters. Lb /Mb also gives us a hint to the star formation 
efficiency in clusters (see below). In any case, the agree- 
ment of brightening between M/Lb and Lb /Mb found 
here can be regarded as indirect support of our conclu- 
sion that single-burst like models seem to be plausible for 
describing the evolution of Lb/ Mb- 

The absolute value of Lb /Mb tells us a hint to the 
star formation efficiency in clusters, i.e., the fraction of 
baryons in clusters used to form stars. If the star forma- 
tion efficiency differs among clusters, the absolute value 
of Lb /Mb would also vary from cluster to cluster. The 
fact that there exist models, such as those with r = 0.1 
Gyr, which reproduce the observed Lb /Mb of many clus- 
ters at different redshifts within the observational errors 
suggest that the star formation efficiency is universal 
among clusters up to z ~ 1. 

In order to see when baryons are converted into stars, 
we plot in figure 7 the predicted evolution of M^,/Mb, 
where the evolution of Mi, is calculated from mass-to- 
luminosity ratios of galaxies predicted by the zp = 3 
models. The values of Mi^/Mb at z = have been nor- 
malized so that they are consistent with the observed 
Lb /Mb of the nearby clusters [|. The predicted values of 
Mi,/Mb at z = are 0.05 — 0.13, depending on r, and are 
consistent with the observed value (the filled square with 
an error bar). This simply implies that the predicted 
(M/Lb)* values at z = fall within (6 ± 3)h{M/LB)Q, 
which is the adopted value of {M/Lb)* for galaxies in 
nearby clusters. As expected, the evolution of Mi,/Mb 
largely differs among the models. A constant M^/Mb is 
predicted by the r = 0.1 Gyr model in the redshift range 
of this figure, while for the r = 5 Gyr model, half of the 
stars present today are formed at z < 1. 

Finally, we mention how to put stronger constraints on 
models using the evolution of Lb /Mb- Unfortunately, the 
steepness of brightening up to z = 1 is not a monotonic 

1 M,/Mt,(z = 0) = (M/LB)r"'(^ = 0) X (LB/Mi,)°b=(2 = o) 
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function of t: the brightening is the steepest for r = 1 — 2 
Gyr models, and models with r = 0.1 Gyr and t — 3 — 5 
Gyr give similar brightening. In order to place further 
constraints using the evolution of L^/Mh, one needs data 
at z > 1: for example, for = 3 and 2, r = 0.1 Gyr 
models predict much steeper brightening at z > 1 than 
r = 5 Gyr models. 

3.2. Comparison with Evolution of Field Galaxies 

The quantity Lb /Mb is the S-band luminosity per 
unit baryon in clusters. The corresponding quantity for 
field galaxies is the ratio of the blue luminosity density 
[Lbq Mpc"*^] to the mean baryon density pb [M^ 
Mpc~^] . In this subsection, we derive Ib / Pb of field galax- 
ies up to z ~ 1 from the literature and compare it with 
Lb/ Mb of clusters. 

Data of Ib are taken (or computed) from recent mea- 
surements of luminosity function based on redshift sur- 
veys: Lilly et al. (1996; CFRS; data points are at 
z = 0.35,0.625,0.875), Ellis et al. (1996; Autofib; z = 
0.085,0.25,0.55), Colless (1998; 2dF; z = 0.11), Loveday 
et al. (1992; APM; z ~ 0.05), Zucca et al. (1997; ESP; 
z ~ 0.1), and Marzke et al. (1998; SSRS2; z ~ 0.025). 
Lilly et al. (1996) give itself while the other papers 
give only luminosity functions in the B band. For those 
except for Lilly et al. (1996), we integrate the luminosity 
function given in each paper from Mb — —25 to —10 to 
obtain ^b- 

In order to compute pb, we adopt VL^ = 0.02/i^'^ follow- 
ing Tytler and his coworkers' results (e.g.. Buries, Tytler 
1998). Their estimates of fi^ are based on measurements 
of deuterium abundance (D/H) of QSO absorption sys- 
tems at high redshifts. Note that the measurements of fif, 
have not completely converged among authors, ranging 
from Vtbh^ ~ 0.01 to 0.02, though Tytler et al.'s measure- 
ments seem to be the most reliable (e.g., Turner 1999). 

The ^B / Pb of field galaxies calculated above are plotted 
in figure 8 as open circles with error bars. A gradual 
increase in Ib / pb with redshift is seen. This is due to the 
brightening oIIb 0- 

A fact which needs attention is that the values of Ib / Pb 
at z < 0.05 are smaller than those at z ~ 0.1 by as large as 
a factor of ~ 2. It is unlikely that the luminosity density 
evolves so rapidly for such a short time from z = 0.1 to 
the present. A possible explanation for this problem is 
that the number density of galaxies in the local (z < 0.1) 
universe happens to be lower than the global value (e.g., 
Marzke et al. 1998), though further investigations are 
needed in order to prove or disprove this explanation. In 

2 For field galaxies, the global star formation rate (p* [A^ yr~^ 
Mpc~^]) has also been measured from UV and emission-line 
luminosities of galaxies (e.g., Madau et al. 1998; for recent 
observations, see Cowie, Songaila, Barger 1999). These mea- 
surements suggest that increases by a factor of ~ 3 — 10 
from z = to 1, but the scatter among the data is still large. 



this paper, we regard the values at z ~ 0.1 as the local 
(z = 0) value. 

3.2.1. Local values for Lb /Mb and Ib/ Pb 

The filled circles in figure 8 indicate the Lb /Mb of the 
distant clusters of galaxies. The filled square at z = 
corresponds to the average of the nearby clusters. We 
find that the local Lb /Mb agrees with the local Ib/ Pb'- 
LB/Mb is (0.024±0.004)(Lb/A^)o and the average of the 
three ?b/p6 values at z ~ 0.1 is (0.026 ± 0.002)(Lb/M)o 

This agreement implies that the blue luminosity per 
unit baryon mass is very close between clusters and fields. 
Then the next question may be whether the stellar mass 
per unit baryon mass (M*/Afb), i.e., the star formation 
efficiency, is the same between clusters and fields (Note 
that morphological type mix largely differs between clus- 
ters and fields and that (M/Lb)* of galaxies varies with 
morphology). In order to examine this, we do a sim- 
ple (but more detailed than that given in §2.1) esti- 
mation of Mi,/Mb of clusters and fields below. We as- 
sume the (M/Lb)* of elliptical and SO galaxies to be 
^h{M/LB)Q and that of spiral and irregular galaxies to 
be 4/i(M/Lb)o 0- Using these values and taking account 
of the mean morphological type mix of the 12 clusters, we 
obtain Mi, /Mb = 0.114 ±0.020. A similar calculation for 
the field galaxies gives A'/^/A/b = 0.097 ± 0.009 (We use 
the type mix given in Colless 1998: /b(E -|- SO)/ZB(tot) = 
34.5% and /b(S + Irr)//B(tot) = 65.5%). Hence, the 
Mi, /Mb of clusters is in agreement with that in fields 
within errors. Renzini (1997) has also obtained a similar 
result from a rough calculation for Mi, /Mb of the Coma 
cluster and field galaxies. 

3.2.2. Evolution of LB/Mb andlB/pb 

From figure 8, it is found that the evolution of LB/Mb 
as a function of redshift is the same as for Ib/ Pb within 
the observational errors. Both 'brighten' by a factor of 

2 — 3 from z = to z = 1. 

The agreement between the evolution of LB/Mb and 
/b/p6 may be interpreted as the mean star formation 
history of cluster galaxies being similar to that of field 
galaxies. However, we think that this agreement is prob- 
ably superficial. The observed global star formation rate 
of field galaxies increases by a factor of ^ 3 — 10 from 

3 A note should be added. The agreement found here does not 
hold if very different values for f2i, and h are adopted, because 
(i) Is/pb changes linearly with and (ii) the dependence of 
Ib/ Ph and L-g/Mf, on h is different {Ib/ Ph h \i f2f,h^ is fixed, 
while Lb/Mi, oc h^'^), though we think that the values of Qi, 
and h adopted in this paper are the most probable at present. 

4 The value for elliptical and SO galaxies is based on van der 
Marel's (1991) study, which found that the M/Lb of elliptical 
galaxies with 1 X 1O1°/i~^Lb0 is 8.4/i(M/Lb)0 . The value for 
the Galactic disk, M/Lb — 3(M/Lb)0, is adopted as the value 
for spiral and irregular galaxies (4/i ~ 3 for h = 0.7). 
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the present epoch to z = 1, and then has a peak at 
z = 1 - 2 (e.g., Madau et al. 1998; Cowie et aL 1999). 
Though this star formation history, which reproduces the 
observed evolution of I Pb as well, should be a solution 
for the mean star formation history of cluster galaxies, 
quite different models can also be solutions, as has been 
seen in §3.1. Simple models having just two parame- 
ters (t and zf) were examined in §3.1, and many models 
have been found to reproduce the observed evolution of 
Lb /Mb, and the color of galaxies in nearby rich clusters 
suggests that r = 0.1 Gyr models are favored (see §3.1). 

In any case, we cannot give a clear conclusion about 
the mean star formation history of cluster galaxies on 
the basis of the current data. As was mentioned in §3.1, 
data at z > 1 are useful to place further constraints on 
the star formation history. More desirable may be data 
of ultraviolet luminosity per unit baryon mass, Ljjy/Mb, 
from which one can measure star formation rate (and 
efficiency) directly. 

4. Conclusions 

We have derived Lb /Mi, for massive (Mgas at the Abell 
radius is > 1 x 10"'^^/i-~^'^M»i) clusters of galaxies up to z ::::: 



We have made a simple estimate of star formation ef- 
ficiency (M^/Mb) to find no difference between clusters 
and fields. To place further constraints on the mean star 
formation history of cluster galaxies needs new data at 
higher redshifts or direct measurements of star formation 
rate. 
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Table 1. Distant Clusters. 



name 


z 




h\ 


n': 


radius'^-' 


error*^ 


A1413 


0.14 


0.028 


7.6 


2.42 


0.77 


21 


A1689 


0.18 


0.042 


14.0 


2.80 


0.78 


21 


A2218 


0.18 


0.043 


8.4 


1.81 


0.42 


33 


A2163 


0.20 


0.037 


3.0 


0.74 


0.26 


48 


CL0500-24 


0.32 


0.037 


3.9 


0.95 


0.53 


30 


CL0939+47 


0.41 


0.089 


1.8 


0.18 


0.14 


— 


RXJ1347.5-1145 


0.45 


0.053 


65.7 


11.3 


1.09 


37 


CL0016+16 


0.55 


0.044 


50.6 


10.4 


1.67 


— 


AXJ2019+1127 


1.01 


0.075 


5.6 


0.68 


0.30 


30 


A2390 


0.23 


0.030 


11.3 


3.3 


0.60 


— 


MS0302+16 


0.42 


0.055 


5.1 


0.84 


0.60 


— 


MS0440+02 


0.20 


0.027 


2.5 


0.84 


0.60 


— 


MS0451+02 


0.20 


0.021 


6.3 


2.7 


0.60 


— 


MS0451-03 


0.54 


0.036 


14.7 


3.7 


0.60 




MS0839+29 


0.19 


0.054 


4.5 


0.74 


0.60 




MS0906+11 


0.17 


0.066 


9.2 


1.3 


0.60 




MS1006+12 


0.26 


0.028 


6.4 


2.0 


0.60 




MS 1008- 12 


0.31 


0.051 


11.2 


2.0 


0.60 




MS1224+20 


0.33 


0.066 


7.1 


0.98 


0.60 




MS1358+62 


0.33 


0.062 


11.4 


1.7 


0.60 




MS1455+22 


0.26 


0.021 


4.9 


2.1 


0.60 




MS1512+36 


0.37 


0.028 


4.2 


1.4 


0.60 





a) In units of /i° '5(Lb/M)o. 

b) In units of X 10"^®. 

c) In units of /j-^-^ x lO^^M®. 

d) Radius in units of Mpc adopted to measure Lb 
and Mgas. 

e) Relative error (%). 
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Figure Captions 

Fig.l. Ly plotted against Mgas for Arnaud et al.'s (1992) 
nearby clusters. The filled and open circles indicate clus- 
ters with and without type mix data, respectively. The 
dashed and solid lines indicate the best fit of a power law 
to all and massive (Mgas > 1 X lO^^/i-^-^M®) clusters, 
respectively. The dotted line corresponds to the best fit 
of L\ oc Mgas to the massive clusters. 

Fig. 2. Ly/Mg^s plotted against the fraction of luminosity 
emitted from elliptical and SO galaxies to the total lumi- 
nosity Lv(E/SO)/Lv for Arnaud et al.'s (1992) clusters 
having type mix data. The filled and open circles indicate 
Mgas > 1 X 10^3/^-2.5^^ ^j^^ M^^^ ^ly. iQi^h-'^-^M^ 

clusters, respectively. 

Fig. 3. L^/Mt as a function of Mgas for Arnaud et al.'s 

(1992) clusters. The filled and open circles indicate clus- 
ters with and without type mix data, respectively. 
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cate the 22 distant clusters and the filled square is for the 

average Ls/Mb of the nearby clusters. The open circles 
with error bars correspond to Is/pb of field galaxies. 
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Fig. 4. Lb versus Mgas for the 12 nearby and the 22 dis- 
tant clusters adopted in this paper. Filled circles indi- 
cate the nearby clusters. Open circles and crosses are 
for 0.1 < z < 0.4 and 0.4 < z < 0.7 clusters, respec- 
tively. The star corresponds to the most distant cluster 
AXJ2019+1127 at 2: = 1.01. The thick solid line indi- 
cates the best fit of Lb Mgas to the nearby clusters. 
The thin solid line and the dotted line correspond to a 
similar fit to the distant clusters at 0.1 < 2; < 0.4 and 
0.4 < z < 0.7, respectively. 

Fig. 5. Dependence of the measurement of LB/Mgas on ilg 
and Ao. The ratio of LB/Mgas(f^o, Ao) to LB/Mgas(0.2, 0) 
is plotted as a function of redshift for (fiojAo) = (1;0) 
and (0.2,0.8) cases. 

Fig. 6. Observed Lq/Mi, of clusters as a function of red- 
shift. The filled circles present the 22 distant clusters 
and the filled square indicates the average Lb/M^ of the 
nearby clusters. Model predictions are overlaid. Pre- 
dicted values of Lb /Mb are normalized to match the ob- 
served value at 2 = 0. Panels (a), (b), and (c) are for 
Zf = 00,3, and 2, respectively. Thick and thin solid 
lines indicate models with r = 0.1 and 1 Gyr, respec- 
tively. Dotted, dashed, and long-dashed lines correspond 
to models with r = 2, 3, and 5 Gyr, respectively. 

Fig. 7. The ratio of stellar mass to baryon mass, plotted as 
a function of redshift. The five lines indicate predictions 
of models with different r. All models are for zp = 3. 
The meanings of lines are the same as in figure 6. The 
filled square corresponds to the observed value. 



Fig. 8. iB/pb of field galaxies as a function of redshift, 
compared with Lb /Mb of clusters. The filled circles indi- 
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